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We investigate the dynamics of football matches. Our goal is to characterize statistically the 
temporal sequence of ball movements in this collective sport game, searching for traits of complex 
behavior. Data were collected over a variety of matches in South American, European and World 
championships throughout 2005 and 2006. We show that the statistics of ball touches presents power- 
law tails and can be described by g-gamma distributions. To explain such behavior we propose a 
model that provides information on the characteristics of football dynamics. Furthermore, we discuss 
the statistics of duration of out-of-play intervals, not directly related to the previous scenario. 



PACS numbers: 02.50.-r. 89.90.+n 01.80.+b, 



I. INTRODUCTION 

The statistical analysis of physical systems has been 
fundamental for the identification of underlying mecha- 
nisms. This kind of analysis has also found applications 
in other fields, such as biology, economics and social sci- 
ences, with a reciprocal feedback for understanding phys- 
ical systems. In any context, when dealing with complex 
systems, where individual, collective and aleatory fea- 
tures may be present, a central interest is to trace general 
statistical properties of the dynamics. This suggests a 
direction to investigate, amongst other human activities, 
collective sports such as the most popular one: football. 
In fact, in football matches, player actions range from 
elementary individual reactions to elaborated strategies 
involving several players, motivating the search for traits 
of complex behavior. 

Amongst the diverse football games, we restrict our 
study to male official association football (soccer). Be- 
cause of its popularity and widespread diffusion in the 
media, an abundant source of observational data is ac- 
cessible. Previous works have dealt mainly with macro- 
scopic features measured over ensembles of matches (cups 
or championships) [l|, H, H, 0], such as the statistics of 
goals. Meanwhile, the present goal is to characterize a 
microscopic dynamics throughout each match. From a 
related perspective, the detection of temporal patterns of 
behavior has been pursued before [5j. Differently, in this 
work, we analyze the stream of ball events throughout 
a match from a statistical point of view. We focus on 
the temporal aspects, without taking into account any 
spatial counterpart, e.g., players or ball trajectories in 
space. 

The results we report in this paper derive from infor- 
mation collected over twenty six matches in South Amer- 
ican and European championships in 2005-2006 (Table 
QJ. Temporal series were obtained from the sequence of 
touches in each football match. 

After the presentation of collected data in Sec. II, we 



expose in Sec. Ill our modeling of the statistics of times 
between touches, starting with simple exponential distri- 
butions, refining the description through gamma distri- 
butions and finally through generalized gamma distribu- 
tions. Basically, we show that a non-stationary Poisson 
process allows to describe with success the main statis- 
tical properties. Before making final observations (Sec. 
V) , we also discuss the statistics of in and out-of-play in- 
tervals (Sec. IV). On one hand our results show the pos- 
sibility of applying statistical physics methods to study 
collective aspects of sports such as football, evidencing 
interesting features. On the other, we exhibit a con- 
crete example where the mechanisms of the "superstatis- 
tics" recently proposed by Beck and Cohen Q, in con- 
nection with Tsallis statistics [7], and where a generalized 
gamma distribution (g-gamma) plays a central role, ap- 
ply. Therefore, our results may provide insights on a 
more general context. 



II. DATA ACQUISITION AND PRELIMINARY 
ANALYSIS 



Data were acquired from TV broadcasted matches, 
with the aid a computer program that records (with pre- 
cision of 10 _2 s) the instants of time at which predeter- 
mined keys are clicked. Thus, time was recorded by click- 
ing when a given player action takes place. Considering 
that the average human reaction time is about 0.1-0.2 s, 
this may lead to a systematic error of that order in data 
recording. However, in the present study we are inter- 
ested in time differences, therefore, any systematic error 
will be reduced. Moreover, some of the matches were 
first recorded and then played in slow motion for data 
collection, showing no significant differences from those 
obtained in real time. Also, in order to check that other 
biases were not being introduced by the observer, the 
authors independently obtained the data for some of the 
matches and verified that the corresponding distributions 
essentially remained the same. 



2 



For each match considered, we monitored the occur- 
rence of touches (kicks, headers, shots, throw-ins, etc.) 
that change the player in possession of the ball. That is, 
touches from one player to himself were not taken into 
account. The instants of occurrence of ball touches were 
recorded without making distinction on the type of touch, 
nor on the subsequent movement of the ball (rolls, flies, 
etc.), nor on whether it was intentional or accidental. 
Moreover, touches were not distinguished by teams. We 
also recorded the time at which each sequence of touches 
ends, that typically corresponds to the instant when the 
whistle is blown. 

First, let us consider the variable T corresponding to 
the time that elapses between two consecutive touches 
occurring without interruption of the match (inter-touch 
time). A typical time series of inter-touch times in a 
match is exhibited in Fig. [TJ This plot manifests the 
discontinuous nature of football activity where the se- 
quences of touches are interrupted by events such as the 
ball leaving the field, player fouls, defective ball, external 
interference or any other reason to stop the game. Then, 
time series are characterized by sequences of ball-in-play 
fragments. 



the exponential decay time r', since in all cases the ini- 
tial decay is exponential. Calculation of r' was performed 
through an exponential fit to the cumulative distributions 
of increments over the interval AT < 4s. The fitting 
procedure, that will denominate WLS along the paper, 
was weighted least squares, with weights w\ = at- 
tributed to each data point (xi,yi), and was performed 
by means of commercial software Origin. Fig. [3Ja) al- 
ready exhibits qualitatively the narrow statistical diver- 
sity found from one match to another, concerning the 
most frequent events. Constraints, such as rules of the 
game and human effort limitations, confine single real- 
izations to a narrow spectrum. In fact, values of r' fall 
within the narrow interval r' = (1.89 ± 0.18) s. Collapse 
of the re-scaled data is observed up to |AT| ~ 2r' only 
(Fig. Elb)). Naturally, deviations from the mean behav- 
ior occur for large AT because events are rare. 



III. MODELING INTER-TOUCH TIME 
STATISTICS 
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Let us consider the random variable X t representing 
the number of occurrences during the period [0,t). If 
the stochastic process {Xt, t > 0} were purely Poissonian 
(as commonly considered for arrival time statistics), with 
expected rate 1/r, then the PDF of inter-touch times 
should be the exponential fr{t) = e~*/ r /T Q. Moreover, 
given two independent variables T 1; T 2 with the same 
exponential distribution, the increment AT = T 2 — T\ 
has the so called double exponential or Laplace PDF 



FIG. 1: Time series of T (time that elapses between 
consecutive touches). The match is Chelsea x Barcelona 
in 22/Feb./2006 (UEFA Champions League). Symbols are 
joined by linear segments whenever there is no interruption 
between touches. 



Typical histograms of inter-touch times T and incre- 
ments AT (between consecutive inter-touch times) are 
presented in Fig. [^a-b). Unless otherwise stated to the 
contrary, in this and subsequent analyses, touches of goal- 
keepers were not considered, because of the singular role 
in the game. These plots, built for the match of Fig. [TJ 
are quite similar to those observed for other matches. In 
all cases the decay of the probability density function 
(PDF) of inter-touch times is approximately exponential 
and the PDF of increments presents a "tent" shape in the 
semi-log representation, suggesting a double exponential 
decay. 

Histograms of the increments, built for several 
matches, are exhibited in Fig. [3Ja). Cumulative distri- 
butions were considered to reduce fluctuations. In order 
to quantitatively characterize each match, we calculated 
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FIG. 2: Normalized histograms of (a) inter-touch times and 
(b) their increments, for the match of Fig. [JJ The inset in 
(a) is a semi-log representation of the same data plotted in 
the main frame. Dashed lines represent: (a) Gamma PDF 
given by Eq. (T2J| with parameters /3 = 3.43 ± 0.17 and r/s= 
0.85 ± 0.05, and (b) the corresponding PDF of increments 
given by Eq. Q. The full line in (b) corresponds to Eq. (|6}. 
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FIG. 3: Cumulative distribution of increments. In (a) 
F\AT\{t) = Prob(\AT\ > t) was computed for each of the 
twenty matches listed in Fig. [4] In (b) , increments are scaled 
by the respective exponential decay time r' obtained from the 
exponential fit to each curve in (a), performed for values up 
to approx. It' only. 



/atW = ^ Wr - (1) 

Although the distribution of increments AT is in good 
accord with a Laplace PDF, at least for central values, 
the distribution of inter-touch times T is clearly not a 
pure exponential (see Fig. (Ha)). Therefore, the expo- 
nential, with one single fitting parameter, constitutes a 
very coarse model for the histograms of inter-touch times. 

Instead, the time interval between touches can be 
thought to be composed, in average, of a certain num- 
ber (3 of independent phases. If each of the phases is 
T-exponentially distributed, then one obtains the Erlang 
distribution: 

Mt) = ^(*/r) /3 - 1 e- t / T = r AT (t), (2) 

defined for t > 0. This PDF is also known as gamma 
distribution, for real f3 > 1. In the particular case (3=1, 
one recovers the pure exponential distribution. However, 
in the vicinity of the origin the PDF of inter-touch times 
has a shape compatible with (3 > 1. In fact, very short 
inter-touch times are not frequent since players are not 
typically so close to each other. Very long inter-touch 
times are also scars since teams dispute ball possession 
almost all the time. Moreover, let us remark that the 
Erlang distribution is commonly used for modeling the 
distribution of times to perform some compound task, 
such as repairing a machine or completing a customer 
service [9]. Also in the present case, when a player re- 
ceives the ball, it is common that he executes more than 



one task, such as, keeping possession of the ball, avoid- 
ing opponents, and passing the ball to another member 
of his team. In what follows we do not restrict (3 to take 
integer values. Then (3 can be interpreted as an average 
number of phases. 

Fig. [2ja) shows the results of WLS fitting the gamma 
PDF to the numerical histograms of inter-touch times. 
We observe a clear improvement in the description of the 
statistics of inter-touch times, in comparison with the 
pure exponential model, for small and moderate times. 
Assuming a gamma distribution, parameters (3 and r 
were determined by means of a least square fit to em- 
pirical PDFs, using statistical weights in order to ponder 
the tails of the distributions. Even so, although a good 
description of numerical histograms is obtained for small 
and moderate values, there are important deviations at 
the tails. This also suggests that the gamma distribution 
may not be a very good model for the present data. 

Fig.[4]exhibits the values of parameters (3 and r for sev- 
eral matches, together with the value of r' (for the cumu- 
lative distribution of increments). Values of r are found 
within the interval [0.68,1.06], with average value 0.91, 
while (3 values fall within the interval [2.6, 4.2], with mean 
3.24. There is a tendency that faster games (smaller r) 
are characterized by a larger 0, as clearly observed in 
Fig. |U There is also a trend that more decisive matches 
or matches played by highly ranked teams have smaller 
r and larger j3, e.g., PA, SL are final matches, AR, CB 
are usual matches of the UEFA champions league. These 
tendencies are expected because in such matches players 
usually save no efforts and strategies are more elaborated. 
If we assume that successive inter-touch times are inde- 
pendent identically gamma-distributed variables, then, 
the PDF of increments becomes 

p-\t\/T POO 

fA T (t)= t20 2 J o [x(\t\+x)f- l e- 2 ^dx. (3) 

Notice that this function has the same asymptotic be- 
havior as the Laplace PDF but it is smoother at the 
origin. In particular, for (3 = 3 (integer value closer to 
the average one) , one gets 

f AT (t) = ^[3 + (3+\t\/r)\t\/r}e-^\ (4) 

Once obtained the values of parameters ((3, r) for the 
distribution of inter-touch times, they were used to pre- 
dict the PDF of increments by means of Eq. ([3]). From 
the results, shown in Fig. [2Kb), one concludes that the 
prediction of the PDF of increments, assuming indepen- 
dence of consecutive inter-touch times, is satisfactory ex- 
cept at the tails. 

Let us discuss some points that may be responsible for 
deviations from the simple Poissonian framework. First, 
we investigated the assumption of independence of occur- 
rences in non-overlapping intervals. We calculated auto- 
correlation functions both for variables T and AT, taking 
into account the intrinsically discontinuous nature of the 
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time series. Hence, only pairs of times belonging to the 
same sequence of touches were considered. No significant 
correlations were detected in the series of T, although the 
series of AT typically presents traits of antipersistence. 
Hence, despite the existence of strategies and patterns in- 
volving several players (some in cooperation and others 
in opposition), the lack of significant correlations indi- 
cates that memory effects in the succession of touches 
are very weak. 

Another possibility for the failure of the simple Poisso- 
nian picture concerns stationarity and homogeneity. In 
order to investigate this aspect, we analyzed for each 
match the number of events N as a function of time, as 
illustrated in Fig. \5\ Beyond the discontinuity of the time 
series, temporal inhomogeneities throughout a match are 
common. First of all, average rates computed over each 
half of the matches are almost always different. In most 
of the cases, these average rates are larger in the first half, 
as soon as players are usually more tired in the second 
half of a match. At a finer time scale, small segments with 
different rate (slope) can be identified (especially during 
the first half, in the case of Fig. [Ha)). This feature is 
a manifestation of the change of rhythms throughout a 
match. We estimated local rates A = 1/r as A = n/t IN , 
where i/jy is the duration of each full sequence of inin- 
terrupted touches and n is the number of touches in that 
sequence. The histogram of rates A is shown in Fig.[5jb). 
Meanwhile panel (c) displays A as a function of time, 
putting into evidence its fluctuating character. 

We will see that the fact that the rate of occurrence A is 
not constant, but instead it is a fluctuating quantity, may 
explain the behavior of the tails of inter-touch time dis- 
tributions. In effect, let us interpret the PDF Tp i/\(t) as 
the conditional PDF fx\\ of variable T given A, where A 
is a stochastic variable. Moreover, let us also assume that 
A is gamma distributed, i.e., f\(x) — T a<K (x). Although 
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FIG. 4: Parameters /3 and r, obtained from the WLS fit 
of the gamma distribution to the histograms of inter-touch 
times, for the matches indicated in the abscissa axis (see Ta- 
ble U for letter code). Parameter r', obtained from cumula- 
tive histograms of increments is also exhibited. Matches were 
ranked by decreasing /3. 
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FIG. 5: Number of events as a function of time (a), for 
the mathc of Fig. [T] Histogram of local rates A (b). Solid 
lines represent WLS fits of a gamma distribution T a ,K with 
(q:,k/s _1 ) = (11.1,0.032) Local rate A as a function of time 
(c). 



this may be only a crude estimation of the distribution of 
local rates, it takes into account its main features, except 
deviations at the tails (see Fig.^b)). Under the assump- 
tions above, the marginal PDF has the form [(| [13, [H[ 



M*)= / dxh{x)f nx {t,x)=m< } - 1 e- t '' 



(5) 



where M is a normalization factor, t = 1/[(q + 

q = 1 + l/(a + /3) and the g-exponential function (e 9 ) for 

negative argument is defined as e~ x = [1 + (q — l)x] 1 -i , 
if q > 1. This PDF that generalizes the gamma dis- 
tribution, is known as F-distribution or also as g-gamma 
distribution [T2|. Panel (a) of Fig.[6]shows the WLS fit of 
the g-gamma function to the same data of Fig. [21 In order 
to asses the goodness of fit we applied the Kolmogorov- 
Smirnoff test [HI]. We calculated confidence levels a by 
determining the largest deviation between the cumula- 
tive distribution that arises from WLS fit and the ob- 
served one. We obtained higher confidence levels for 
the g-gamma model. As illustration, for match CB, the 
gamma and g-gamma fits yielded a = 7% and 32%, re- 
spectively. Furthermore, the chi-square value of fit and 
the correlation coefficient for match CB were (x 2 ,i? 2 )= 
(0.0009, 0.99) (against (0.002, 0.97) for the simple gamma 
distribution) . 

Although at the cost of introducing one more parame- 
ter, the g-gamma model is satisfactory for the full range 
of values. The advantage of introducing one more param- 
eter was quantified through Akaike information criterium 
(AIC=2fc — 21nL, where k is the number of parame- 
ters and L the maximum likelihood) and also through 
Schwarz criterium (SIC=/cln(n) — 21nL, where n is the 
number of observations), which penalizes more strongly 
the introduction of free parameters [l4| • In ah cases the 
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FIG. 6: Normalized histograms of inter-touch times, (a) 
for the match of Fig. [Tj and (b) for all the matches inves- 
tigated. Dashed lines correspond to the gamma distribu- 
tion with parameters (J3,t/s)= (3.43 ± 0.17, 0.85 ± 0.05) (a) 
and (3.12 ± 0.12, 0.95 ± 0.05) (b). Solid lines correspond 
to the g-gamma distribution with parameters (/3, r/s, q)= 
(4.45 ± 0.35,0.43 ± 0.07,1.066 ± 0.007) (a), (4.20 ± 0.17, 
0.43 ± 0.04,1.076 ± 0.003) (b). Insets: linear representations 
of the same data. 



g-gamma distribution yielded lower values. For exam- 
ple, (AIC, SIC)= (3726, 3741) (against (3760, 3772) ob- 
tained with the simple gamma) for CB and (82694, 82718) 
(against (83650, 83666)) for the global set. 

Furthermore the PDF of increments, which generalizes 
Eq. ©, namely, 

p oo 

/atW« [x(\t\+x)f- 1 e-^ T e-^+ x ^ r dx , (6) 
Jo 

also describes better numerical data (see Fig. [3(b), where 
the full line is the predicted PDF of increments using the 
parameters of Fig. [6(a)). 

In Fig. [5(b), all the matches were merged. In this case, 
the merging procedure by itself might give rise to the dis- 
persion of A responsible for the behavior of the tails. This 
global analysis, however, is useful for characterizing the 
average behavior of football activity as a whole. Since, 
as observed before, diversity is not very high amongst 
matches, one observes for the merged data a behavior 
similar to that observed for a single match (illustrated in 
Fig. [5(a)), although the statistics is poorer in the later 
case. 

The introduction of parameter q allows to describe bet- 
ter the statistics of rare events. Notice that, in compari- 
son with the gamma fits, the g-gamma ones yield (3 about 
one unit larger and r about one half smaller. Alterna- 
tively, assuming that the fluctuating rates obey a gamma 
distribution r QjK , one can obtain q (i.e., q = l + l/(a+(3)) 
and t (i.e., r = l/[(a + for the resulting g-gamma 

distribution. The values of q are very close to those 
obtained by directly fitting the g-gamma distributions 
(q ~ 1.07). Whereas, the resulting values of r are larger 
than those obtained from g-gamma fits but still of the 
order of Is. Therefore our model is selfconsistent. The 



lack of a complete matching of parameter r is due to 
diverse reasons. On one hand, there is a certain degree 
of arbitrariness in the definition of local rates, that in 
our case were computed over each continuous sequence 
of touches, through n/tiN. Moreover, the distribution 
of local rates as here denned are not strictly gamma, 
but only approximately. Finally, also parameter /3 is an 
averaged quantity since the number of phases may fluc- 
tuate throughout a match. Nevertheless, the comparison 
between theoretical and empirical distributions supports 
the present model as a better approximation than the 
simple gamma distribution. 

IV. IN AND OUT-OF-PLAY INTERVALS 

In Fig. [7(a) we present the histogram of tjN (dura- 
tion of the intervals without interruption). Notice that 
tjN = 52i=i^i- Then, its PDF can be obtained as 
p{tiN) — J2 n >i ptfiN \n)P(n), where P(n) is the distri- 
bution of the number of touches n in each continuous 
sequence. The conditional PDF is in first approxima- 
tion r^ jnr , assuming that the n inter-touch times are in- 
dependent identically r^ iT stochastic variables. On the 
other hand, P(n) follows approximately the exponential 
law e- n ' n °/n (Fig. [7(b)), being n a ~ 7.54 ± 0.21 (from 
WLS fit), while (n) ~ 8.4 ± 0.3. Under the assumptions 
above, one obtains a PDF that can be well approximated 
by an exponential with characteristic time tin — (n)f3r, 
consistent with the numerical results in Fig. [7(a). 
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FIG. 7: Statistics of in-play intervals. Cumulative histogram 
of tiN (a). The straight line corresponds to a exponential 
with characteristic time tin = 28s. Normalized histogram of 
the number of touches in each continuous sequence (b). 

The histogram of times elapsed between sequences (du- 
ration of intervals when the ball is out-of-play) was also 
computed and it is exhibited in Fig. [5(a). Since both tiN 
and touT statistics are poor over a single match (about 
one hundred events), we merged the records of several 
matches (AV,RA,MB,VZ,LF,BM,CB, following Table H}. 
In these cases goal-keeper movements were computed. 

Although the statistics of out-of-play intervals appears 
to display a power-law behavior, a more careful analysis 
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points to a log- normal statistics (Fig. [SJb)). Up to now, 
the results were well understood within the framework 
of non-homogeneous Poisson processes. However, in the 
case of out-of-play intervals, the PDF basically obeys a 
log-normal statistics that is not straightforwardly related 
to the previous scenario. 



V. FINAL OBSERVATIONS 

The statistics of touches can be understood on the ba- 
sis of Poissonian arrival or point processes. Basically, 
events are not simply Poissonian but can be thought 
as composed of different phases, as also observed for 
compound tasks. This explains the behavior of the his- 
togram of inter-touch times in the vicinity of the ori- 
gin. Meanwhile, non-homogeneities in the rate of occur- 
rence (changes of rhythm) throughout a match appear 
to be responsible for the power law tail in the distri- 
bution of inter-touch times, giving rise to a g-gamma 
function. It is noteworthy that a similar mechanism 
based on compound distributions for the obtention of 
such PDFs has been proposed within the context of the 
"superstatistics" [gjj, where the fluctuating parameter is 
the temperature. Here, we provide an example where 
g-gamma distributions arise as a consequence of the fluc- 
tuating nature of a relevant parameter. 

All the main features here exhibited and discussed for 
the matches in Table Q] arc also observed, in a prelim- 
inary analysis, for the sixty four matches of the 2006 
World Cup (results not shown). Within the simple Pois- 
sonian framework, the effective characteristic time r (of 
the order of one second) does not change significantly 
from one match to another. Meanwhile, parameter (3 ex- 
hibits a greater variation, remaining approximately be- 
tween 2 and 4, with average number of tasks close to 
3. In general there is a tendency that r is shorter and 
(3 larger in more decisive matches or matches played by 
highly ranked teams. These trends are qualitatively ex- 
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FIG. 8: Statistics of out-of-play intervals. Cumulative his- 
togram of tour (a). In (b), normalized histogram of In tour 
(symbols) and Gaussian fit with mean 2.85 and variance 0.45 
(solid line). 

pected as far as in such matches players usually save no 
efforts, playing faster and using more elaborated strate- 
gies. In fact, larger (3 already suggests more developed 
or complex actions. The introduction of a further pa- 
rameter, q, which reflects the degree of inhomogeneity 
of rhythms, improves the description of the tails. The 
statistics of the duration of sequences of touches, inter- 
rupted by fouls, ball leaving the field, etc., can also be 
derived within the same approach used for the statistics 
of touches. 



On the other hand, the statistics of intervals between 
sequences of touches is of a different nature, belonging to 
the log-normal class. There are diverse mechanisms that 
may give rise to such statistics. As an example, for time 
series observed in turbulent flows, it has been attributed 
to a multiplicative random process fill ]. However this 
issue should be further investigated and deserves separate 
work. 
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Match(cup) 



BR Barcelona lxl R.Madrid(S) 

RB R.Sociedad 0x2 Barcelona(S) 

AV A.Bilbao lxl Villareal(S) 

RA R.Madrid 3x0 Alaves(S) 

VB Valencia 1x0 Barcelona(S) 

MB Mallorca 0x3 Barcelona(S) 

VZ Valencia 2x2 Zaragoza(S) 

RM Rcggina 1x4 Milan (I) 

LF Leverkusen 2x1 E. Frankfurt(G) 

BM Borussia M. 1x3 B. Munchen(G) 

CB Chelsea 2x 1 Barcelona (E) 

AR Arsenal 0x0 R. Madrid (E) 

SL Sao Paulo 1x0 Liverpool(C) 



Date 



02/Apr./06 
19/Mar./06 
26/Feb./06 
19/Feb./06 
12/Feb./06 
29/Jan./06 
19/Jan./06 
12/Feb./06 
28/Jan./06 
27/Jan./06 
22/Feb./06 
08/Mar./06 
18/Dec./05 



Match(cup) 



PA Sao Paulo 4x0 Atletico-PR(L) 

DC D. Cali 0x1 Corinthians(L) 

CU Corinthians 2x2 U. Catolica(L) 

PD Palmciras 2x0 D. Tachira(L) 

PB Parana 2x0 Botafogo(B) 

PC Parana 1x2 Corinthians(B) 

CS Corinthians 1x3 Sao Paulo(B) 

SA Santos 2x0 Atletico-PR(B) 

SF Sao Paulo 1x0 Flamengo(B) 

PS P. Santista 0x5 Sao Paulo(P) 

SP Santos 1x0 Palmeiras (P) 

SC S. Caetano 2x1 Corinthians(P) 

SM Santos 3x2 Marilia(P) 



Date 



14/Jul./05 
15/Feb./06 
22/Feb./06 
25/Jan./06 
03/Aug./05 
14/Mar./06 
07/May/06 
23/Apr./06 
16/Apr./06 
12/Feb./06 
05/Mar./06 
08/Feb./06 
22/Jan./06 



TABLE I: Matches of South- American and European championships. Cups are B (Brazilian) , C (FIFA World Club) , E (UEFA 
Champions League), G (German), I (Italian), L (Libertadores), P (Brazilian, Sao Paulo State), S (Spanish). 
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